Modeling electrostatic patch effects in Casimir force measurements 
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Electrostatic patch potentials give rise to forces between neutral conductors at distances in the 
micrometer range and must be accounted for in the analysis of Casimir force experiments. In 
this paper we develop a quasi-local model for describing random potentials on metallic surfaces. 
In contrast to some previously published results, we find that patches may provide a significant 
contribution to the measured signal, and thus may be a more important systematic effect than was 
previously anticipated. Additionally, patches may render the experimental data at distances below 
1 micrometer compatible with theoretical predictions based on the Drude model. 

PACS numbers: 31.30.jh, 12.20.-m, 42.50.Ct, 78.20.Ci 



I. INTRODUCTION 

The Casimir effect QUI] is a remarkable consequence 
of vacuum field fluctuations which, in its simplest man- 
ifestation, leads to the attraction of two neutral ideal 
conducting plates. At very short distances quantum fluc- 
tuation forces dominate the interaction between neutral 
objects making them an essential consideration for micro- 
electro mechanical devices (MEMS) and atom traps, 
among others. The comparison between experimental 
measurements and theory for Casimir forces between 
metallic plates has been a matter of debate in recent 
years. This debate is of particular importance if this 
comparison is used to derive constraints on hypothetical 
new short-range interactions appearing in addition to the 
gravity force in unification models [SHE]. 

Two recent experiments are at the heart of this debate. 
Casimir force measurements by the IUPUI group [51 IIP), 
performed at distances smaller than 750 nm, were inter- 
preted by the authors as excluding the dissipative Drude 
model and agreeing with the lossless plasma model. This 
has led to a discrepancy between experiments and phys- 
ically motivated theoretical models, such as the Drude 
model, for real conductors which exhibit dissipation. In 
distinction, a recent experiment by the Yale |llj was able 
to measure Casimir forces at distances up to 7 /xm and 
was interpreted by the authors as being in agreement 
with the Drude prediction, including quantum as well as 
thermal fluctuations, once an electrostatic patch contri- 
bution has been taken into account. 

It is known that patch effects are a source of concern 
for Casimir experiments [T2HT6] , as well as for other pre- 
cision measurements [I7H26J. For the Yale experiment 
the patches were assumed to be much larger than the 
gap D between the spherical and planar plates used in 
the measurement. Under these conditions the patch force 
is found to be proportional to RV t 2 ms /D in the proximity 
force approximation (see below), where R is the radius 
of curvature of the spherical plate and V Tms is the root- 
mean-square (rms) voltage of electrostatic patch poten- 
tials [TT]. For the IUPUI experiment, a patch analysis 



was performed with different assumptions leading to the 
conclusion that the patch effect had a negligible influ- 
ence Unfortunately, it was not possible in any of 
these experiments to measure the patches independently. 
It follows that the conclusions of the theory-experiment 
comparisons heavily rely on the patch models used in the 
data analysis. 

In this paper, we revisit electrostatic patch effects and 
analyze their possible influence in Casimir force measure- 
ments. Our approach is based on the method pioneered 
by Speake and Trenkel [T2] with the electrostatic patches 
described in terms of a power spectral density. However, 
we will develop a model for the power spectral density 
differing from the one proposed in |12| and used in [5I|10|. 

Our model is based on the observation that bare metal- 
lic surfaces are composed of crystallites, each of which 
constitutes a single patch, where the local surface volt- 
age is determined by the local work function [37j . By as- 
suming that through the surface preparation process the 
crystallographic orientation, and hence the correspond- 
ing work function, of each crystallite is determined inde- 
pendently and randomly we can infer that voltage corre- 
lations are restricted to points lying on the same patch: 
We refer to this as quasi-local correlation. Our model 
with quasi-local correlations can be compared to the case 
of quenched charge disorder in dielectrics [3U], and 
also shows close similarities with models proposed re- 
cently to describe patch correlation functions for atomic 
or ionic traps [3TJ [32] • 

We will show that the voltage correlation function from 
our quasi-local model strongly differs from that initially 
proposed in [12] and used in [10] . As a result, in con- 
trast to the claims of [HI [TO], patches may have a sig- 
nificant contribution to the IUPUI measurements. In 
addition we will qualitatively address the issue of sur- 
face contamination which is expected to lead to larger 
correlation lengths and reduced voltage fluctuations [25J . 
Given that the degree of contamination is unknown we 
perform a fit of the patch model we propose to the differ- 
ence between measurements and the Casimir theoretical 
prediction based on the Drude model and find that it 
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qualitatively explains the residual signal. For the Yale 
experiment, our results will essentially reproduce those 
obtained in [TT] . 



II. ELECTROSTATIC PATCH EFFECT 

In the present section, we recall a few general results 
of interest, assuming that the validity conditions of the 
proximity force approximation (PFA) are satisfied, that 
is, the radius of the sphere used in the experiment is much 
greater than the sphere-plane distance. In this case, the 
expression for the force gradient G sp (derivative with dis- 
tance of the force F sp ) in the sphere-plane geometry is 
written as follows in terms of the pressure P pp (the force 
per unit area) calculated between two planes 



G sp {D) 



dF sp {D) 
dD 



2irRP pp (D). 



(1) 



This expression is used throughout the paper for both 
Casimir and patch effects. 

The basic description of the patch effect after [12] is a 
statistical ensemble of patch potentials Vi(r) on the sur- 
faces of two planar plates labeled i — 1,2. The potentials 
are assumed to have zero mean (Vi(r)) — 0, and to be de- 
scribed by the two-point potential correlation functions 



C«(r) = (Vi(r)^-(O)) = 



4^' 



r dj [k] 



(2) 



In the plane-plane geometry, points on the planes are 
denoted in cartesian coordinates as r = (x, y) and the 
point is an arbitrary origin. The correlation functions 
Cy(r) and therefore the power spectra Cy[k] are also 
assumed to be isotropic. The relations between these 
two functions can be written 

Cij(r) = I dkk J (kr) C tJ [k], 
27T Jo 

Cij[k] = 2tt / drr J a (kr) C y (r), (3) 
Jo 

where we have simply denoted r = |r| and k = |k| and 
where J n (x) is the n-th order Bessel function [33 . The 
patch power spectrum CV,- [k] corresponds to the notation 
Cij(k) in |12j . As usual, the variances and covariances 
are given by the integrals 



Ctf(O) = (ViVj) 



1 f°° 

— / dkkCij[k] 
2t Jo 



(4) 



The pressure due to electrostatic patches in the plane- 
plane geometry can be computed exactly [12] as 



ppatc hp) = 



4tt 



dk k 3 
sinh 2 (kD) 



(5) 



x {Cu[k] + C 22 [k] - 2C 12 [k] cosh(kD)} . 

It is worth emphasizing at this point that the integral 
is reduced to a very simple expression when patch sizes, 



with a typical value denoted £ pa tch, are larger than the 
distance D. In this case, all wavevectors k contributing 
to the integral ([5| satisfy kD <C 1, so that the pressure 
scales universally as 1/ D 2 , irrespective of the particular 
details of the power spectrum (Eq. Q is used) 
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{C 11 [k]+C 22 [k}-2C 12 [k}} 

D < ^patch- (6) 



The above result is expected from the analogy with a 
parallel plate capacitor with prescribed voltages. In con- 
trast, when the relevant wavevectors no longer satisfy 
the above inequality, different models for the patch power 
spectrum result in different predictions for the patch con- 
tribution to the pressure. 

It is also worth mentioning here some conditions for 
the expression §5§ of the electrostatic patch pressure be- 
tween two plates to be valid. A fundamental assumption 
in this analysis is that the ergodic hypothesis is satisfied, 
which means that the distribution of patches within the 
interaction area is a fair approximation of the ensemble- 
averaged distribution function defined by the power spec- 
trum C-ij [k] . When applied to two plane plates of finite 
area A, we expect this assumption to be well satisfied if 
the effective interaction area contains a large number of 
patch correlation areas A 3> ^p atc h- F° r t ne sphere-plane 
geometry, the effective area of interaction is of the order 
of nDR, leading to the validity requirement 
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patch' 



(7) 



In the following two subsections we recall a model used 
in [H] and [HI HO], and introduce another model with 
quasi-local correlations which we think to be a better 
description of sputtered surfaces. 



A. The sharp-cutoff model 

We now discuss the model of patch correlations which 
was proposed as an example in |12j and then used in 
[9l [10] to assess the contribution of electrostatic patches 
to the Casimir force measurements. 

It is a simple description based upon two assumptions: 
a) the power spectrum of patches is an annulus in k- 
space possessing no other dependence than a sharp cutoff 
at small (fc m i n ) and large (A; max ) wavevectors (hence the 
name sharp-cutoff model); b) there are no cross correla- 
tions between the two plates {C\ 2 = 0) . This model gives 
the power spectrum for a single plate as 



AttV 2 

Cu[k] = 



k 2 

"'max 



-0(fcmax-fc)0(fc-fcmin), (8) 



where V 2 ms is the variance of the potential on one plate 
and 9 is the Heaviside step function. 

In order to determine the parameters of this model, the 
authors of [S] used the further assumptions : c) based on 
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AFM images of the surfaces, the minimum and maxi- 
mum grain sizes of the samples were determined to be 

^"ch = 25 nm and ^atah = 300 nm ! d ) the P atch sizes 
were assumed to be the same as the grain sizes and the 

cutoffs in fc-space were derived from the inverse maximum 
and minimum grain sizes A: m j n = lixji™^^ — 20.9/xm -1 
and fcmax = 27r/£™ 1 " ch = 251/im" 1 ; e) the rms voltage 
was obtained by computing the variance of the work func- 
tions over the different crystallographic planes of gold, 
which led to V^ ms ~ 80.8 mV. Using the five assump- 
tions a) to e), it was concluded in [9] that the patch 
pressure had a negligible influence on the estimation of 
the Casimir force. A reasonable agreement was then ob- 
tained between the experimental data and the prediction 
for the Casimir pressure using the lossless plasma model 
(more discussions below). 

Now we will argue that model ([8]) is not a good de- 
scription for the patch power spectrum for the surfaces 
used in the experiments, and later on, we will also ques- 
tion the relation between patch and grain sizes. In order 
to make the former point clear, let us write the correla- 
tion function Ca (r) of patches in real space which can be 
obtained through an inverse Fourier transform ^ from 
the spectrum ^ 



C«(r) = 2V£ 



^max^l (^max^") ^min^l (^minT) 



(k 



2 

max 



(9) 



As one moves away from coincidence the correlation func- 
tion Ca (r) oscillates between positive and negative values 
with a period of the order of the smallest patch size, and 
is contained within an envelope decaying as r~ 3 / 2 (see 
FigjTJ. These oscillations imply that the patch potential 
shows correlations as well as anti- correlations in space. 
Such behavior could be expected for surfaces exhibiting 
some kind of antiferroelectric ordering (where the con- 
figurational energy is minimized when adjacent surface 
dipoles are antiparallel), but will unlikely describe the 
random potentials on sputtered surfaces. 

As already stated, the strict relation between patch 
sizes and grain sizes, assumed in the analysis of [S], has 
also to be questioned. The adsorption of contaminants 
on the surfaces alters patch sizes which, as a result, do 
not necessarily correspond to the grain sizes [28] , We 
expect that contamination leads to an effective smearing 
of the patch layout, so that patch sizes will be larger than 
grain sizes while the voltage variance will be less than the 
value obtained for a clean sample from the assumption 
e) discussed above. 



B. The quasi-local correlation model 

We now propose another patch model which we think 
to be a better motivated description of the patch corre- 
lation function for the surfaces used in the experiments. 

To model the layout of cystallites on a plate, we choose 
a random patch layout and afterward assign a random 
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FIG. 1. Comparison of the sharp-cutoff and quasi-local patch 
models described in subsections II-A (dashed lines) and II-B 
(solid lines), respectively. Plot (a) shows the voltage correla- 
tion functions in real space while plot (b) shows the associated 
spectrum in Fourier space. All plots correspond to the corre- 
lation function C = Ca divided by Kms- On the lower plot, 
the sharp-cutoff spectrum discussed in II-A is multiplied by 
a factor of 100 in order for it to appear at the scales shown. 
The parameters used for both models, discussed in subsec- 
tions II-A and II-B, are taken from |5j, but however, do not 
correspond to the same average patch size. 



potential to each patch. For a given micro-realization of 
patches we write the voltage over the whole surface as 



V(x) 



u a e a (x). 



(10) 



The sum is over all patches, v a is a random variable de- 
scribing the voltage on patch a, and the function o (x) 
is defined to be 1 for x on the ath patch, and otherwise. 

We now obtain the two-point voltage correlation func- 
tion by performing ensemble averages over all micro- 
realizations of the patch voltages and layouts. Physically, 
the voltage on each site is determined by the crystallite 
face exposed to the surface. As we assume that each 
crystallite is deposited with a random crystallographic 
orientation and that each deposition is statistically inde- 
pendent we can infer that 



{v a Vb) 



S ab V; 



2 

ms 5 



(11) 



where the expectation value (...)„ averages over the volt- 
age fluctuations only and S a b is the Kronecker delta. 
Also, note that we are implicitly assuming that there 
are no cross correlations between the patches on differ- 
ent plates, C12 = 0. Using (10 1 and (111 we construct 
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the two-point voltage correlation for a single micro- 
realization of the patch layout 

(y(x)v(x')) u = 7 r L £ e Q (x)e a (x'). (12) 

a 

The final step in constructing the ensemble-averaged 
voltage correlation function is to average over all patch 
layouts. We carry this out by exploiting several symme- 
tries: 

1. We assume that the patches are distributed uni- 
formly and isotropically which implies that the av- 
erage patch associated with any given point on the 
surface is circular with a radius determined from a 
distribution of patch sizes. In reality no patch is 
circular and this notion of patch radius should only 
be taken in a statistical sense. 

2. For any two points on the sample surface, the 
voltage correlation function C(x, x') is propor- 
tional to the number of patches which contain both 
points (among all micro-realizations). By employ- 
ing the statistical description of patches, as de- 
scribed above in 1, the correlation will be computed 
by summing over all circular patch centers and sizes 
as depicted in Fig. [2] 

3. As a check one can verify that the correlation 
at coincidence is the constant V 2 ms . Moreover, 
translational and rotational invariance implies that 
C(x, x') depends only on r = |x — x'|. 




FIG. 2. The voltage correlation function, C(x, x'), is con- 
structed by summing over all circular patches which contain 
both x and x'. This is undertaken by integrating over patch 
centers y, and accounting for the distribution in patch sizes 
with the distribution II(^). 



Given these considerations we find the following form 
for the correlation function: 

/>oo 

C(x,x')= / diU(£) (13) 
Jo 

J d\j 9(1/2 - |x - y\)8(e/2 - \x' - y|), 



W 2 



where the integral over y, constrained by the ^-functions, 
sums over all patches of size I which contain both points. 
The final integral over i averages over patch sizes where 
II(£) is the distribution of patch diameters. Note in par- 
ticular that the translational invariance of the correlation 
function is made apparent by the change of variables 
x' — y — > z . Subsequently performing the integration 
over y reveals the rotational invariance of the final result 



C(r) = C u (r)-- 



2V 2 



di n(f) 



cos 



1 



(14) 



The patch power spectrum can then be obtained 
through the Fourier transform pi). Some interesting 
properties can be given at this point. First, the integra- 
tion of the correlation function over all space is simply 



C[k = 0] = 7 Tr£ 2 V x 



(15) 



with £ 2 the variance of the distribution H(£). Second, the 
integral of C[k] over wavevectors is just the variance of 
the potential J Q °° dkkC[k] = 2nV^ ns . 

One can derive some universal scaling laws for the 
patch contribution to the pressure in some limiting cases. 
When the patches are much larger than the gap (D <C ^) , 
the expression Q is obtained. This 1 / D 2 scaling law for 
the pressure (and the corresponding l/D for the energy 
per unit area for planar plates) is universal for all patch 
power spectral densities whenever the typical patch sizes 
are much larger than the gap. In particular, this scaling 
was used in [TT] to model the patch effect [34]. In the 
opposite limit, where the typical patch sizes are much 
smaller than the gap, one can obtain a simple scaling law. 
In this case the spectrum C[k] is approximately constant 
over the wavevector range k ~ l/D which provides the 
most significant contribution to §5^. We then find, when 
using ( 15 ), 
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(16) 
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We emphasize at this point that this 1 / D 4 scaling law is 
generic for all spectra having a finite limit at k = 0, but 
does not hold when C[k] vanishes at k = 0. In particular, 
in the model discussed in II-A, there is a sharp-cutoff of 
the power spectral density at fc m i n > 0. In this case, the 
pressure |5]) is exponentially small when k m \ n D 3> 1, that 
is also D ^> ^™^h- The leading order contribution indeed 
comes from the exponential tail of 1/ sinh 2 (£;£>) and is 
much smaller than the result found in the generic case 
( 16 1. This point will play a crucial role in the comparison 



to experimental data discussed in the next section. 

Before entering this discussion we choose a specific 
form for the patch size distribution n(^) which is similar 
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in spirit to the sharp-cutoff modei discussed in subsection 
If-A. By assuming the patch sizes are distributed uni- 
formly within a finite interval between a minimum ^™" ch 
and maximum ^J^* h value, the probability distribution 
is 



U(£) = 



Max tfmin 

patch patch 



and has the following moments 



Max 
patch 



/?mm 
patch 



(17) 



(18) 



max \2 
patch/ 
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min \2 
patch / 



/7m ax pram 
patch patch 



Additionally, we would like to remark that we have also 
considered other size patch size distributions H(£) (log- 
normal, Gaussian, generalized gamma, etc.) and have 
found similar results for the pressure in all cases. 

We emphasize that, despite some similarity in the con- 
struction of the two models discussed in subsections II- A 
and II-B, they correspond to very different correlation 
properties, the most striking difference resulting from a 
nonvanishing value for C[k = 0] in the quasi- local model 
which gives a distinct large distance behavior. In par- 
ticular, a patch model employing quasi-local correlations 
was recently adopted to describe heating in ion traps and 
dissipation in cantilevers [31 j . There, the observed large 
distance (D ^> J) scaling of electric field noise (oc D~ A ) 
is linked with a nonvanishing value of C[0]. 

To estimate the effects of contamination, we will as- 
sume that the patch power spectrum on a dirty surface 
takes the same form as on a clean surface (i.e., also given 
by the quasi- local model) , with the exception that the pa- 
rameters of the model are altered by the contaminants. 
Let us stress that quasi-local correlations may be not as 
accurate for contaminated surfaces as for clean ones. We 
employ the above assumptions in a preliminary manner 
to account for the properties of contaminated surfaces, to 
be confirmed by dedicated studies to come in the future. 



III. COMPARISON WITH EXPERIMENTS 

We now compare the theory and experiments by calcu- 
lating the Casimir force from the Drude model, and the 
patch pressure arising from the model with quasi-local 
correlations. To make the comparison we first calculate 
the plane-plane Casimir pressure P pp (D) at temperature 
T using the Lifshitz formula [55H57] . We use tabulated 
optical data for gold , extrapolated to low frequencies 
with a Drude model to describe the contribution of con- 
duction electrons, e CO nd(w) = 1 — flp/(uj(iu + 17)), where 
flp is the plasma frequency and 7 quantifies the damping 
rate. To account for roughness corrections to the Casimir 
pressure we adopt the simplest formulation based on an 
additive scheme (Eq. (33) in [9]). We will call the re- 
sulting pressure as the "Drude model" Casimir pressure 
(D). 



As already stated, we use the PFA to relate the experi- 
mental data corresponding to the sphere-plane geometry 
to the predictions calculated in the plane-plane geome- 
try, for the Casimir and the patch effects. In the IUPUI 
experiments the sphere-plane force gradient G sp is mea- 
sured, which is related to the equivalent plane-plane pres- 
sure as in (JlJ. In the Yale experiments the sphere-plane 
force F sp is measured, which is related similarly to the 
plane-plane energy per unit area. 

After subtracting from the experimental data the the- 
oretical predictions for the Casimir interaction, we find a 
residual signal 



6P Dmdc (D) ee P^ pcrimcnt (£>) - P° rudo (D). 



(19) 



The question we address in the following is whether or 
not the residual §p Urudc can be explained by a reason- 
able modeling of patch effects. The criterium is then to 
minimize the remaining difference between the residual 
signal and the patch pressure 6P Bludc (D) - pP atch (D). 
The residual is defined here for the Drude model and 
may be as well be defined for the plasma model. The 
patch pressure pP atch (D) is then defined for a given patch 
model, say in particular the sharp-cutoff (subsection II- 
A) or quasi-local (subsection II-B) models. 



A. Data analysis for the IUPUI experiment 

For the comparison with the IUPUI experiment we 
compute the Casimir force at room temperature T = 295 
K using tabulated optical data extrapolated to low fre- 
quencies with a Drude model with parameters flp = 8.9 
eV for the plasma frequency and 7 = 0.0357 eV for the 
damping rate. Root mean square roughness heights for 
the plane and the sphere are 3.6 nm and 1.9 nm, respec- 
tively. These permittivity and roughness parameters are 
the ones reported in |10j . 

We collect in Fig(3] the information needed to com- 
pare IUPUI experimental data with predictions from the 
Drude model and modelings of the patch effect. We plot 
the residuals SP Bludc defined as in ( 19 ) as points with 



pDrudc 
PP 



error bars and the patch pressure P pa tch for different 
patch models as lines. The error bars represent the to- 
tal experimental error described in Fig. 2 of [TU] at 67% 
confidence. The theoretical predictions for the Casimir 
pressure Pp p ludc are calculated for the Drude model as 
described above and assumed to have no error. There 
are four different patch models represented in Fig|3] : 

1. The solid curve is the estimation of the patch ef- 
fect using all the assumptions of subsection II-A. 
The patches are thus described by the sharp-cutoff 
model ^ with the parameters fc max = 251 ^m _1 , 
fcmin = 20.9 /im _1 and Vi. ms = 80.8 mV (these are 
the parameters used in [5]). 

2. The dotted curve is the result of the quasi-local 



correlation model ( 14 1 with the patch size distri- 



6 



bution (17) described in subsection II-B. The pa- 
rameters, ^p a '" ch = 25nm, ^tdi = 300nm and 
Kms = 80.8 mV, correspond to the assumptions 
that the patch sizes are given by the grain sizes 
and the rms voltage is determined by the variance 
of the work function over the different crystallo- 
graphic planes (these are the same parameters used 
in item 1 above). 

3. The long-dashed curve is obtained from a 
least-squares minimization of the difference 

jpDrudc ( D ) _ ppatchp^ usmg the quas i_l oca l 

patch correlation model given by (14) and (17 1. 



As ^™ a '" ch is found to have a small influence, we 
fix it to the smallest grain size ^™ a " ch = 25 nm as 
discussed above. The best fit on the two remaining 
parameters gives ^™ a a t * h ~ 2476 nm and V rms « 9.2 
mV and results in qualitative agreement between 
the residual and the fitted patch pressure. The 
associated first moments of the patch size distri- 
bution are I = 1251 nm and P = (1437 nm) 2 . The 
reduced-x 2 for this fit, calculated using the total 
error bars from Fig 2. of [TO] at 67%, is 0.814. 
It is important to note that the values of the fit 
parameters and quality of the fit are very sensitive 
to the sample's optical parameters, in particular 
to the plasma frequency used in the extrapolation 
of optical data to low frequencies [39]. However, 
one should avoid giving too much importance to 
any of these values of reduced-x 2 as a measure 
with statistical significance of experiment-theory 
agreement. Indeed, the influence of sample de- 
pendency of optical parameters, the use of a very 
crude description of roughness corrections to the 
Casimir pressure, and, most importantly, the lack 
of precise information of the patch correlation 
function in actual experimental samples, all imply 
that the fits obtained with the quasi-local model 
for patches have a qualitative nature; dedicated 
patch effects measurements are required to make 
metrological claims (sec the Conclusions for further 
discussions). 

4. The short-dashed curve (underneath the long- 
dashed curve) is a fit of a phenomenological model 
proposed by Carter and Martin [32]. The correla- 
tion function of this model, based on a Monte-Carlo 
simulation of patch layouts, can be expressed in 
terms of a shifted Gaussian and is specified by the 
rms voltage and the average patch area w 2 , related 
to our patch radius via w ^fit 1/2. Our best fit 
values are £ w 1229 nm and Vr ms = 8.6 mV with 
reduced-x 2 of 0.812. 

After this description of the information gathered on 
Fig. [3| let us now comment on the significance of the 
various results: 

1. The solid curve reproduces and confirms the cal- 
culations which were performed to quantify patch 
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FIG. 3. Comparison of the residual <5p Drudc between the ex- 
perimental pressure in [TO] and the Drude prediction (points 
with error bars at 67% confidence taken from Fig. 2 of [10] ) 
with patch pressure P pa tch for four different patch models 
(more details in the main text) : 1. The solid curve is the 
result of the sharp-cutoff model (with assumptions of subsec- 
tion II-A) ; 2. The dotted curve corresponds to the quasi- 
local patch correlation model assuming that the patch sizes 
are given by the grain sizes and that the rms voltage is given 
by the variance of the work function over different crystallo- 
graphic planes ; 3. The long-dashed curve is the result of a 
best-fit on the parameters (^ max and F rms ) of the quasi-local 
patch correlation model ; 4. The short-dashed curve (under- 
neath the long-dashed curve) is a fit of a phenomenological 
model proposed in [32] . The inset shows the residual sig- 
nal resulting from subtracting the fit of the quasi-local model 
(long-dashed curve) from 5p Drudc . 



effects in [HI HO]- With the assumptions described in 
subsection II-A, the calculated patch pressure is in- 
deed far too small to explain the difference between 
experimental data and theoretical predictions using 
the Drude model. 

The dotted curve gives the result of the quasi-local 
model of patch correlations ( [17] ) with parameters 
determined as was done in |9l IIP), but here for a 
different patch spectrum model. As a striking il- 
lustration of the importance of this difference, the 
calculated patch pressure is now larger than the dif- 
ference between experimental data and theoretical 
predictions using the Drude model. This illustrates 
the highly model dependent nature of the computed 
patch pressure. Thus, patches may be an impor- 
tant systematic effect for which their contribution 
to the measured signal should ideally be assessed 
independently of any Casimir force measurement. 

The long-dashed curve corresponds to a least 
squares fit of the quasi-local correlation model to 
the residual SP Bludc . With the best- fit parameters 
ppatch an d Kms, this model qualitatively fits the dif- 
ference between experimental data and theoretical 
predictions using the Drude model. These param- 
eters have reasonable values: ^atch ^ s l ar S er than 
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the maximum grain size on the samples, and V Ims 
smaller than the rms voltage for a clean sample 
[HI QJj] . This suggests the presence of contaminants 
on the sample surfaces [25] . 

4. The best-fit of the phenomenological model pro- 
posed in [32] is essentially indistinguishable from 
that of the quasi-local correlation model (long- 
dashed curve). The best-fit values for 1 and V ims 
are consistent with the average patch size and rms 
voltage obtained from the best-fit parameters of the 
quasi-local model. 

At this point, we also want to comment on the validity 
requirement ([7]), which allows one to calculate the patch 
effect in the sphere-plane geometry within the PFA. This 
requirement ensures that the effective area of interaction 
between the sphere and the plane, of the order of irRD 
for a sphere of radius R, contains a large number of el- 
ementary patch areas, so that the sum over the micro- 
realization of patches on a given plate is a good effective 
description of the statistical ensemble-average given by 
the power spectral density. With the numbers in |10j . 
that is a radius of curvature of the sphere R — 151. 3^m 
and a shortest distance D min — 160nm, the interaction 
area is ttRD w 76(/im) 2 . Meanwhile, the average patch 
area is (ir/4:)£ 2 m 1.6(/im) 2 , there is a large number of 
elementary patch areas (sa 48) within the effective area 
of interaction, but it is possible that one could expect a 
small correction to the patch pressure at short distances 
when the ergodic hypothesis begins to break down. 

For completeness we have also studied the residual 
SP piasma (D), as defined in Eq. (TTob , with the excep- 



tion that we have compute the plane-plane Casimir pres- 
sure PPjJ asma (Z)) using the "plasma model", instead of 
the Drude model. More precisely, we have computed the 
pressure using for the permittivity e(i£) the "generalized 
plasma model" : 



fi 2 

: - plasma (ze) = i + 



6 

■E 

3 



fj 



(20) 



where the first two terms correspond to the permittivity 
for the plasma model for conduction electrons (dissipa- 
tion of conduction electrons is set to zero ad hoc with- 
out physical justification), and the second sum of terms 
accounts for the interband transitions of gold [ID]. To 
account for roughness corrections to the Casimir pres- 
sure we use the same additive scheme employed above. 
Computing 8P plasma (D) in this way, we have confirmed 
the findings of [TU], namely that a negligible contri- 
bution of the patch effect leads to an agreement of data 
with theoretical predictions using the plasma model. We 
note, however, that the patch pressure calculated from 
the quasi-local model, with sizes and voltages used in 
[9], is much larger than the difference between the mea- 
surements and the plasma prediction, as shown in Fig. 
[4] [J2]. We think that this result constitutes a serious 



\ 




200 



300 



400 500 

D(nm) 



600 



700 



FIG. 4. Comparison of the residual ,5pP lasma (points with 
error bars at 67% confidence taken from Fig. 2 of [10]) with 
patch pressures given by the sharp-cutoff model and the quasi- 
local model: 1. The solid curve is the result of the sharp-cutoff 
model (with assumptions of subsection II- A). We find, con- 
sistently with the analysis in [9] [10] , that the patch pressure 
from this model gives a negligible contribution to the mea- 
sured signal. 2. The dotted curve corresponds to the quasi- 
local patch correlation model adopting the same parameters 
used in [5] [TO] for the patch size and the rms voltage. In 
distinction to the sharp-cutoff model, we find that the quasi- 
local model gives a large signal as compared to the residual 

^pplasma 



warning against the claims according to which the plasma 
model would be confirmed with a high confidence level 
by Casimir experiments performed with real metals [44] . 



B. Data analysis for the Yale experiment 

In addition to analyzing the IUPUI experiment we now 
apply the same models to the recent experiment by the 
Yale group [TTJ. A patch analysis was already carried 
out in [TT] and it led to a good agreement between ex- 
perimental data and the Drude model. This analysis only 
considered the asymptotic form oc 1 /D of the plane-plane 
energy due to patches Here we extend the analysis 
by using the more general expression ([5| for the patch 
pressure with the quasi-local patch correlation function 
described in subsection II-B. Because we have no infor- 
mation regarding grain or patch sizes in the Yale experi- 
ment, we will focus our attention on best-fit estimations 
of the parameters ^™ a a * h , ^™ c h; an d Kms characterizing 
the quasi- local patch correlation function ( 14|17 ). 

To analyze Yale experimental data we first compute 
the Casimir force using tabulated optical data extrapo- 
lated to low frequencies with the Drude model using the 
plasma frequency Qp = 7.54 eV and the dissipation rate 
7 = 0.052 eV employed in [TT]. We set the temperature 
to be T = 295K. The roughness correction to the Casimir 
force is ignored as it gives a negligible correction to the 
force at the distances considered in the Yale experiment. 
Fig(5] shows the difference of the Yale experimental force 
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FIG. 5. Comparison of the residual SF Dludc between the 
data of QT] and the computed sphere-plane force with the 
associated patch pressure i?P atch . The dots correspond to 



SF 



with the errors bars including only the experimental 



error in the force determination. The solid line is a best-fit of 
the patch force within the quasi-local model of subsection II- 
B. The inset shows the corresponding residual ( 5J , P la8ma with 
the same convention employed in the main figure. 



data and the Casimir force prediction us ing the Drude 
model, SF Drude (defined by analogy with (19)), depicted 



by points with error bars (we assume no error for the 
theory). The solid curve shows the resulting patch force 
for parameters arising from a least-squares minimization 
of the quantity <5i^ Drudo — ^ppatch us j n g ^he quasi-local 



patch correlation model (17 1 described in subsection II- 
B. The best fit parameters are given by ^™ a ^h ~ 614/zm, 
^paTch ~ 566/im, (corresponding with £ = 590^m) and 
Vrms ~ 3.9 mV. We should point out, however, that the 
result of the best-fit is essentially insensitive to the details 
of the patch power spectrum. Indeed, since the residual 
^pDrude m ^ ne Yale experiment has an approximate l/D 
power law we can infer using ^ that the typical patch 
size is much larger than D for the whole range of dis- 
tances explored in the experiment (0.7 /im - 7/im). Per- 
forming a constrained fitting by requiring that ^™ a ^* h be 
less than some predetermined value (e.g. 500 ^m), yet 
still satisfying the constraint £ ^> D, we were able to 
verify that a good fit can still be achieved over a large 
range of patch sizes. In summary, we point out the result 
of our fitting using the more detailed quasi-local patch 
model confirms the patch treatment in [TT] . 

Finally, we also report for the sake of completeness 
some supplementary test we performed for comparing the 
data in [TT] with the predictions of the plasma model de- 
scribed by Eq.(20) (see the inset of Fig. [5]). The good 
agreement obtained for the Drude model is dramatically 
degraded. Therefore, we confirm the result obtained in 
[llj that patches cannot explain the difference between 
the experimental data and the plasma model in Yale 
data. 



IV. CONCLUDING REMARKS 

In this paper, we have analyzed the patch contribution 
to Casimir experiments with a model featuring quasi- 
local voltage correlations. Our model is derived from 
well- motivated physical principles and shares key fea- 
tures with experimentally verified patch models used to 
describe ion trap heating and cantilever damping [31j . 
Thus, for the description of the surfaces used in the ex- 
periments discussed in this paper, we believe that this 
model is more appropriate than the sharp-cutoff model 
which has been used to the same aim in previous publi- 
cations [91 [TO]. 

Due to the large difference in the patch power spec- 
trum, in particular for small wavevectors, the quasi-local 
model gives a larger contribution than the sharp-cutoff 
model. As a striking consequence, when the patch sizes 
are deduced from the grain sizes (as was done in [91 [TO]), 
the quasi-local model produces a patch pressure larger 
than the difference between the experimental data and 
the Drude (and plasma) Casimir prediction, whereas the 
sharp-cutoff model produces a negligible patch pressure. 
Therefore, it is important to emphasize that because 
of the combination of: a) the highly model dependent 
nature of the computed patch pressure and b) the po- 
tentially large patch contribution to the measured sig- 
nal, patches may lead to nonnegligible systematic ef- 
fects. This necessitates an independent measurement of 
patch effects in order to meet metrological standards for 
Casimir force measurements. 

We have also used the new quasi-local patch model to 
fit the difference between experimental data of the IUPUI 
experiment [5] HUJ and the theoretical prediction for the 
Casimir pressure. The latter was computed taking into 
a) tabulated optical data extrapolated to low frequencies 
by means of the Drude model, and b) roughness effects 
modeled by a simple additive technique. We have found 
best-fit parameters for the average patch size and for the 
rms voltage that are consistent with a contamination of 
the metallic surfaces, which is expected to enlarge the 
patch sizes (with respect to grain sizes) and smear the 
patch voltage (with respect to those of a surface of bare 
crystallites) [28) . Indeed, surface contamination is ex- 
pected, and we believe that preferential adsorption [35] 
and saturation of contaminants may be compatible with 
the observation of reproducible results in experiments re- 
peated several times with different samples [43] , 

Taken together, our results constitute a strong warning 
against the previously published claims of an agreement 
of Casimir experiments with the plasma model, and an 
elimination of the Drude model 44 . However, we want 
to emphasize that they do not constitute yet a proof of 
agreement of experimental data with the new model. The 
parameters of the patch model have been fitted and it is 
still possible that the qualitative agreement thus obtained 
is a fortunate output of the fitting procedure rather than 
an explanation of the experimental data. 

In this paper we have focused our attention on only the 
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IUPUI and Yale experiments, but of course the analysis 
can be repeated for other Casimir measurements between 
metallic plates as well [ISHSi] . 

A better characterization of the surfaces used in the 
experiments is now key to reaching firmer conclusions. 
The patch distributions can be measured with appropri- 
ate technologies such as Kelvin probe force microscopy 
which can achieve the necessary size and voltage resolu- 
tions [S51[5S]. In addition, the study of cold atoms and 
cold ions trapped in the vicinity of metallic surfaces |22j 
or the role of patch effects in other precision measure- 
ments [2~4"H2"Tj] are other ways for accessing information of 
interest for our problem. Let us repeat at this point that 
our new quasi-local model is similar to recent proposals 
for patch physics used to achieve a better understanding 
of atomic and ionic traps [3U [22] • 

The challenges of forthcoming studies may be stated 
as follows. First, it is important to confirm the hypothe- 
sis that the patch voltages show quasi-local correlations, 
and to better specify the power spectrum which quan- 
titatively describes these correlations. Second, it would 
also be interesting to study how the patch power spec- 
trum depends on contamination, in particular, fabrica- 
tion, treatment, history of the samples, and on tempera- 
ture. Finally, an independent determination of the patch 
power spectrum could lead either to a confirmation of the 
best-fit analysis presented in this paper or to new ques- 
tions. This study is important not only for the test of 
the Casimir effect, a central prediction of quantum field 



theory, but also for the searches of the hypothetical new 
short-range forces predicted by unification models [SHE]. 
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